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Abstract Oscillating solid bodies have frequently been used for studying the 
properties of normal and superfluid helium. In particular, the transition from lam- 
inar flow to turbulence has attracted much interest recently. The purpose of this 
note is to review several central features of this transition in oscillatory flow, which 
have been inaccurately formulated in some recent work. 
PACS numbers: 47.27.Cn, 47.37.+q, 67.25.bf, 67.25.dg 



1 Introduction 

Various types of oscillating bodies have been used for investigating the properties 
of liquid helium: vibrating wires, oscillating disks, grids, spheres and, recently, 
quartz tuning forks [ 1 ] . 

At low velocity amplitudes of the oscillating body the flow is laminar. In this 
regime the amplitude grows linearly with the driving force. In the superfluid the 
drag force is given by the interaction with the quasiparticles and one can distin- 
guish between a hydrodynamic and, at low temperatures, a ballistic interaction, 
depending on the ratio of the quasiparticle mean free path to the size of the body. 
In case of a sphere, the hydrodynamic drag force is known quantitatively from 
Stokes' formula. For vibrating wires an approximation by a straight cylinder is 
made to allow for a quantitative analysis using Stokes' formula for a cylinder. For 
the more complicated geometry of a grid or a tuning fork the analysis is more 
quantitative. Nevertheless, the tuning fork and the vibrating wire have proven to 
be very useful secondary thermometers because the damping force scales with the 
quasiparticle density. In the ballistic regime the drag is given simply by the scat- 
tering of quasiparticles off some geometrical cross section. Finally, between the 
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hydrodynamic and the ballistic regime slip effects become important and a quan- 
titative analysis is complicated. In the normal liquid the drag force in the laminar 
regime is classical and therefore of lesser interest here. 

At large velocity amplitudes the flow becomes turbulent on either side of the 
oscillating body extending into a volume determined at least by the size and the 
amplitude of the body or even filling the entire measuring cell. In this regime 
the amplitude is found to grow with the square root of the driving force, both in 
the normal liquid and in the superfluid. This signals a drag force that increases 
with the square of the velocity amplitude, similar to classical fully developed tur- 
bulence for steady motion of the body. For a sphere the turbulent drag force in 
superfluid ^He was studied in detail [2, 3]. Turbulent drag originates from the 
superfluid component. Only near the lambda temperature nonlinear drag of the 
normal component begins to contribute but full turbulence in the normal phase 
was not observed. With the tuning forks, however, full turbulent drag in the nor- 
mal liquid was reached at sufficiently high velocities fllH]. The transition from 
laminar oscillating flow to turbulence is of great interest, not only in the superfluid 
regime, where shedding of quantized vortices will occur at some critical velocity 
10], but also in the normal liquid because little is known experimentally about the 
different drag force when compared with steady motion. In the following Section 
2, elementary properties of the classical case will be summarized and compared 
with recent work. This will then be the basis of a discussion of results obtained in 
liquid '^He in Section 3. 

2 Transition to Turbulence of Classical Oscillatory Flow 

There are three independent length scales in oscillatory flow: the size L of the 
oscillating body, the oscillation amplitude A, and the viscous penetration depth 
5 — \j2v jiX), where V is the kinematic viscosity of the liquid and to is the oscilla- 
tion frequency. From these 3 quantities one can obtain 2 independent dimension- 
less numbers. Because the transition to turbulence is due to the nonlinear term of 
the Navier-Stokes equation when the velocity v becomes large, the following 2 
numbers are relevant [7J : 

1. The ratio of the nonlinear term (vV)v ^ v^/L to the viscous term vV^v ^ 
vy /L}. This ratio is the Reynolds number Re = L\/v. Inserting v = (OA gives 
Re^2LA/5^. 

2. The ratio of the nonlinear term to the time derivative v ^ Aco^. This is called 
the Strouhal number Sr in Ref.f/'] which is given here hy Sr =A/L. Note that 
the Strouhal number does not depend on the viscosity and hence remains valid 
in an ideal liquid. 

Thus, the situation is more complicated than in steady flow where only the Reynolds 
number is relevant. Now, one has to take into account the Strouhal number as well. 
There are 2 limiting cases, though, where the situation is simpler 10]: 

1. If L ^ 5, i.e., at low frequencies or low temperatures in superfluid ^He (but 
still in the hydrodynamic regime), where 5 of the normal component diverges, 
the time dependence of v can be neglected. The situation is analogous to sta- 
tionary flow and consequently the nonlinear term is negligible if /?e ^ 1 . 
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2. li 5, i.e., at large frequencies or in "^He close to or above 7\ , the nonlinear 
term can be neglected only if the i'r <C 1 . The viscous term is much smaller 
than Aft)^ and the Reynolds number 2AL/5^ need not to be small. 

It is clear that no pair of numbers other than Re and Sr is relevant for the 
transition to turbulence of oscillatory flow. For example, the choice to consider 
the Keulegan-Carpenter number (Kc) and the Stokes number (St) instead [8], is 
not useful. While Kc = 2kA/L is essentially the same as the Strouhal numbeiQ, 
the Stokes number is St = (Ol? /v = 2L^/5^. Because L/8 may be either a large 
number or a small one in the laminar regime, no information on the transition to 
turbulence can be inferred from its value. The reason for this fact is that St is given 
by the ratio of the Reynolds number 2AL/5^ to the Strouhal number A/L. Hence, 
it does not depend on the nonlinear term, which drops out. It just compares the 2 
linear terms of the Navier-Stokes equation: it is given by the ratio of v AaP' to 
the viscous term vV^v ^ vAft)/L^. 

Moreover, it is also clear that one cannot replace L by 5 in Re and Sr as sug- 
gested recently when discussing experiments with tuning forks on the transition 
to turbulence in liquid helium 14), |5D, see below. This would leave only A and 5 
as independent length scales from which only one dimensionless number, namely 
A/5, can be obtained. This would give Re = 2A/8 and Sr=A/5 which obviously 
does not make sense because by definition both numbers refer to independent ra- 
tios of the nonlinear term to the linear ones, and also Sr would then depend on the 
viscosity. 

3 Transition to TXirbulence in Liquid Helium 

Recent experiments with tuning forks on the transition from laminar flow to fully 
developed turbulence cover a temperature range from 1.3 K up to 4.2 K for the 
liquid and include also gaseous helium at 78 K 1 4, 5]. In the normal liquid and in 
the gas a gradual transition from the linear behavior of v as a function of the drive 
F to the square root dependence is observed, extending over more than 2 orders 
of magnitude in v. This is similar to what is known from steady motion, see, e.g., 
Ref. 1 7]. By extrapolating both regimes a "critical" velocity \cr can be defined 
where the laminar drag force i^j™ = Av becomes equal to the turbulent drag force 

^ira? = yv^, hence \cr = A/y. The drag coefficients k and y = Cdpo /2 can be 
obtained from the data {Cd is a numerical drag factor that depends on the geometry 
of the body, p is the density of the liquid, and a is the area normal to the flow)o 



This "critical" velocity is found to scale as \Jva) |4, 5]. For a sphere the drag 
coefficient A is known analytically (Stokes' solution) and therefore the validity of 
this scaling can be proven |^. For the geometry of a fork, A obviously has a similar 
dependence on V and ft) and differs only by a numerical factor. It is important to 
note, that this scaling behavior is understood without postulating, as was done 

' In classical fluid dynamics often a different definition of the Strouhal number is used and 
the definition here is called Keulegan-Carpenter number. 

^ To be precise, for nonlinear drag forces the principle of energy balance (energy gain from 
drive = energy loss from drag) has to be applied for inferring the drag force from the driving 
force. But this will introduce only a numerical factor of order one in y which is not important 
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in Refs. 10, 0], that the length L in the Reynolds number should be replaced by 
the penetration depth 5. Finally, in view of the wide interval between the laminar 
regime and full turbulence, extending over more than 2 orders of magnitude, it 
appears questionable that this "critical" velocity has any physical meaning. At 
least for steady flow around a sphere or a cylinder it has none that I know of. 

In the superfluid phase of '*He the situation is quite different. Well below 
where the normal fluid density is small compared to the superfluid density, the 
transition to turbulence is sharp, even discontinuous or hysteretic iB [3l lol IToll . 
The reason is that turbulence in the superfluid phase occurs abruptly at a critical 
velocity when vorticity is created. The classical Strouhal number is obviously not 
applicable for the transition to turbulence in a quantum fluid. 

Experiments with a sphere 1 , 1 1] show that the turbulent drag force is no longer 
proportional to but instead has a constant shift F'"/^^ = y - Fq. This implies 
that superfluid turbulent drag sets in at a temperature independent critical velocity 
Vcr = yFo/y- Up to now, no theory exists that could describe the measured critical 
velocities quantitatively, probably because both geometry and surface properties 
of the oscillating body will have an influence. Furthermore, careful experiments 
with a vibrating wire show that the critical velocity may be affected by remanent 
vortices ll^ [l2[l . The normal phase may remain in the laminar regime or gradually 
start deviating from it at high velocity amplitudes. With the sphere, e.g., very 
little nonlinear normal drag is found near 10, [H], which is reasonable because 
Sr < 0.5 in this case. The implication is that quantum turbulence in the superfluid 
phase of "^He does not necessarily lead to classical turbulence in the normal phase, 
in contrast to what seems to be generally believed. 

Experiments with an oscillating grid [13] indicate that the resonance frequency 
is lower when turbulence is produced. This is attributed to an enhanced effective 
mass of the grid due to vortices attached to it. However, with wires the opposite 
behavior is found ITqI . Very recently, a decrease of the resonance frequency 
of the wire is reported^in a vortex free environment, and an increase if remanent 
vorticity is present |12]. (In the latter case the spring constant of the oscillator 
is believed to be increased by vortices attached to the body and the walls of the 
cell.) Skepticism is expressed with regard to the interpretation of the enhanced 
mass effect Isl fTill . With the sphere no change of the resonance frequency due to 
turbulence production was detectable within a resolution of better than 10^^. This 
situation clearly deserves further investigation. 

In the ballistic regime, at very low temperatures, the flow pattern switches 
intermittently from laminar flow (better: potential flow) to turbulence and back 
iSIiilEll]. The lifetimes of the turbulent phases are exponentially distributed and 
can be attributed to statistical fluctuations of the vorticity fl^. Turbulence appears 
to be stable when sufficient power from the drive is available so that the turbulent 
lifetimes exceed the measuring timeEl The lifetimes of the laminar phases, how- 
ever, which have a Rayleigh distribution flT'], remain to be explained. Finally, 
metastable potential flow is observed above y^r whose lifetime is limited only by 



^ Because the mean turbulent lifetime is found to increase exponentially with the power in- 
put from the drive but never really diverges, superfluid turbulence is inherently unstable. It only 
appears to be stable because of the finite observation time of the turbulent phases (up to 100 min- 
utes in Ref. [11]). Interestingly, also in classical turbulence an inherent instability has recently 
been observed in pipe flow H7i1 . 
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natural radioactivity, which may create local vorticity in the superfluid that triggers 
the decay of the potential flow fllll . Also this observation of metastable potential 
flow awaits a theoretical interpretation. 

4 Summary 

Several problems with the interpretation of experiments on the transition to turbu- 
lence of oscillatory flow of ^He have been discussed. Some have been resolved by 
reviewing the classical case. In the superfluid, some remain open for further inves- 
tigation. Corresponding experiments in superfluid ^He remain beyond the scope 
of the considerations here, because the situation is quite different: pair breaking 
and Andreev scattering of quasiparticles by the superflow contribute to the drag 
on the oscillating body, and the strength of mutual friction determines whether or 
not turbulence can exist. 
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